Mon. Not. R. Astron. Soc. 000, [TH?] (2008) Printed 31 October 2008 (MN M£X style file v2.2) 



WASP- 14b: 7.3 Mj transiting planet in an eccentric orbit 

Y. C. Joshi 1 *, D. Pollacco 1 , A. Collier Cameron 2 , I. Skillen 3 , E. Simpson 1 , I. Steele 6 , 
R. A. Street 9 , H. C. Stempels 2 , D. J. Christian 1 , L. Hebb 2 , F. Bouchy 14 , N. P. Gibson 1 , 
G. Hebrard 14 , F. P. Keenan 1 , B. Loeillet 14 ' 16 , J. Meaburn 10 , C. Moutou 16 , B. Smalley 5 , 
3 ;I. Todd 1 , R. G. West 7 , D. R. Anderson 5 , S. Bentley 5 , B. Enoch 4 , C. A. Haswell 4 , 

C. Hellier 5 , K. Home 2 , J. Irwin 12 , T. A. Lister 12 , I. McDonald 5 , P. Maxted 5 , 
q ; M. Mayor 13 , A. J. Norton 4 , N. Parley 4 , C. Perrier 15 , F. Pont 11 , D. Queloz 13 , 
^ :R. Ryans 1 , A. M. S. Smith 2 , S. Udry 13 , P. J. Wheatley 8 , D. M. Wilson 5 



oo 
O 
O 
(N 



1 Astrophysics Research Centre, School of Mathematics & Physics, Queen's University, University Road, Belfast, BT7 INN, UK 
2 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife KY16 9SS, UK 



, 3 Isaac Newton Group of Telescopes, Apartado de Correos 321, E-38700 Santa Cruz de la Palma, Tenerife, Spain 

O ii 4 Department of Physics and Astronomy, The Open University, Milton Keynes, MK7 6AA, UK 

I ' 5 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 

' 6 Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead, CH41 1LD, UK 

<4_> 7 Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH, UK 

22 ' 8 Department of Physics, University of Warwick, Coventry CV4 7AL, UK 

, 9 Las Cumbres Observatory, 6740 Cortona Dr. Suite 102, Santa Barbara, CA 93117, USA 

10 Dept of Physics and Astronomy, University of Manchester, UK 

^v^j 11 Physikalisches Institu, University of Bern, Sidlerstrass 5, 3012 Bern, Switzerland 

• 12 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street MS- 16, Cambridge, MA 02138-1516, USA 

qq ' 13 Observatoire de Geneve, Universite de Geneve, 51 Ch. des Maillettes, 1290 Sauverny, Switzerland 

' 14 Institut d'Astrophysique de Paris, UMR 7095 CNRS, Universite Pierre & Marie Curie, 9S ftls bvd. Arago, 75014 Paris, France 

■^-J- , 15 Laboratoire d'AstrOphysique, Observatoire de Grenoble, Univ. J. Fourier - BP 53, F-38041 Grenoble Cedex 9, France 
16 Laboratoire d'Astrophysique de Marseille, CNRS (UMR 6110), BP8, 13376 Marseille Cedex 12, France 



O 

OO 
O 



Accepted 2008 October 31. Received 2008 October 29; in original form 2008 September 10 



ABSTRACT 

We report the discovery of a 7.3 Mj exoplanet WASP-14b, one of the most massive 
^ . transiting exoplanets observed to date. The planet orbits the tenth-magnitude F5V 

star USNO-B1 11118-0262485 with a period of 2.243752 days and orbital eccentricity 
e = 0.09. A simultaneous fit of the transit light curve and radial velocity measurements 
yields a planetary mass of 7.3±0.5 Mj and a radius of 1.28±0.08 Rj. This leads to a 
mean density of about 4.6 gem -3 making it densest transiting exoplanets yet found 
at an orbital period less than 3 days. We estimate this system to be at a distance of 
160 ± 20 pc. Spectral analysis of the host star reveals a temperature of 6475 ± 100 K, 
log g = 4.07 cms -2 and vs'mi = 4.9± 1.0 kms -1 , and also a high lithium abundance, 
logTV(Li) = 2.84±0.05. The stellar density, effective temperature and rotation rate 
suggest an age for the system of about 0.5-1.0 Gyr. 

Key words: planetary systems: individual: WASP-14b — stars: individual: GSC 
01482-00882 — techniques: photometric — techniques: radial velocity 



1 INTRODUCTION 

The giant exoplanets that transit across the disks of their 
host stars are of great interest due to their impact on our un- 
derstanding of planetary structure. Since the first discovery 

* E-mail: y.joshi@qub.ac.uk 



of a t ransiting exoplanet HD 209458b l|Charbonneau et al. I 
l200d : iHenrv et al.T feoOOV more than 40 transiting systems 
have been found around nearby stars. Transit light curves, 
along with their radial velocity motions, provide a wealth 
of information about the system including precise mass, 
radius and mean density of the planet. This in turn al- 
lows us to probe their internal structure by comparing their 
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Figure 1. The Super WASP discovery light curve for the WASP- 
14 (1SWASPJ143306.35+215340.9). The data points are phased 
using the ephemeris, T (HJD) = 2454463.57583, and folded on 
the orbital period of 2.243752 days. 



physical p arameters with mod e ls of planetary struct ure and 
evolution (jCuirlot et. al Il200fj : iFortnev et al. Il2007n . Given 
the importance of these systems, several wide-field sur- 
veys are in p r ogres s to fin d transiting exoplanets, e .g. HAT 
' Bakos et al l l2004h- XO l|McCullough et al l 120051), TrES 



dO'Donovan et al.N2006r i, and SuperWASP (|Pollacco et all 
l2006l l. 

In this paper, we report the 14th exoplanet discovered 
in the SuperWASP survey which is orbiting around a 10th- 
mag F5-type star in the Northern Hemisphere. The paper 
is organised as follows. In §2, we give details of the discov- 
ery light curve and outcome of our follow-up photometric 
and spectroscopic observations. Spectral analysis to deter- 
mine system parameters is described in §3. The stellar and 
planetary evolutionary status is discussed in §4. Finally, our 
results are briefly summarized in §5. 



2 OBSERVATIONS 

2.1 SuperWASP Observations 

WASP-14 (= 1SWASPJ143306.35+215340.9 = USNO-B1 
1118-0262485 = GSC 01482-00882) is an F5V star with V 
— 9.75 and B — V = 0.46 in the constellation Bootes. It 
was monitored by the SuperWASP instrument on La Palma 
during 2004 May to 2007 June. A detailed description of 
the SuperWASP instrument, observing proc e dure and data- 
processing issues is given in iPollacco et al.l ((2006J), while a 
discussion of the candidate selection procedure is explained 
in Collier Cameron et al. (2006, 2007). A total of 7338 data 
points were obtained for this star in four observing runs and 
in two different cameras. The SuperWASP discovery light 
curve of WASP-14b, as shown in Figured] reveals a transit 
recurring with a period of 2.243752 days, and a total dura- 
tion of about 180 minutes between first and fourth contact. 



2.2 Photometric follow-up 

Initial followup photometric observations of WASP-14 were 
carried out by Las Cumbres Observatory, S. Arizona's 81- 
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Figure 2. (a) Phase folded transit light curves from the Su- 
perWASP, Tenagra and Keele data, (b) A high precision transit 
light curve of WASP-14 observed with the RISE. The best- fitting 
model curves determined from a simultaneous MCMC fit (see de- 
tail in §3.3) are shown by the solid lines. 



cm Tenagra robotic telescope and Keele University Obser- 
vatory 60-cm telescope in I and R bands respectively. All 
the images were bias subtracted and flat fielded before do- 
ing aperture photometry around bright stars. Light curves 
along with the SuperWASP discovery data around the tran- 
sit phase are shown together with model fits in Figure [2j a). 

In order to obtain a higher precision photometr}Q of 
WASP-14, we observed the target field using the high-speed 
imaging ca mera RISE mounted on the robotic 2-m Liverpool 
Telescope (jSteele et al. 1120081 ; Icibson et al. Il2008h . We ob- 
served the transit on the night of 2008 May 04 in single broad 
band V+R filter as part of the Canarian Observatories In- 
ternational Time Programme. The RISE Camera consists of 
a thermoelectrically cooled 1024x1024 pixel CCD detector 
which has a pixel scale of ~ 0.55 arcsec/pixel and a total field 
of view of ~ 9.4x9.4 arcmin. During this period we obtained 
7200 frames of 2 sec exposures in 2x2 spatial binning. Aper- 
ture photometry was performed using the IRAF DAOPHOT 
routine around the target and 3 bright, non- variable compar- 
ison stars using an aperture of radius 8 pixels. We binned 
the images in 30-second blocks to increase signal-to-noise 
ratio. To derive the differential magnitude of WASP-14, we 
produced a ratio of the combined flux from the compari- 
son stars to that from WASP-14. We modelled the transit 
light curve using a Markov-chain Monte-Carlo (MCMC) al- 



The data is available upon request. 
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Table 1. Radial velocities for WASP-14 from FIES/NOT and 
SOPHIE/OHP data. The first six data points were taken from 
FIES instrument while the others were obtained with SOPHIE. 



Time 
(BJD-245,0000) 


RadVel 
(kms -1 ) 


&RV 
(kms^ 1 ) 


Bis Span 
(km s" 1 )" 


4461.7400 


-5.8705 


0.0080 


0.0667 


4462.7710 


-4.1180 


0.0054 


0.0629 


4465.7490 


-4.8489 


0.0078 


0.0332 


4466.7850 


-5.4037 


0.0044 


0.0388 


4490.7700 


-5.6209 


0.0087 


0.0392 


4490.7820 


-5.6383 


0.0081 


-0.0069 


4508.5938 


-5.3563 


0.0112 


0.0410 


4509.5238 


-5.0654 


0.0114 


0.0070 


4509.5830 


-4.8313 


0.0098 


0.0270 


4510.5074 


-4.5874 


0.0129 


0.0540 


4510.5599 


-4.7205 


0.0106 


0.0330 


4510.6509 


-4.8798 


0.0099 


0.0260 


4510.6569 


-4.8954 


0.0107 


0.0110 


4510.6630 


-4.9138 


0.0111 


0.0670 


4510.6860 


-4.9758 


0.0100 


0.0130 


4510.7215 


-5.0906 


0.0108 


0.0320 


4510.7274 


-5.1086 


0.0112 


0.0390 


4511.5572 


-5.6865 


0.0107 


0.0200 


4511.6822 


-5.3425 


0.0114 


0.0370 


4512.5406 


-4.2173 


0.0109 


0.0220 


4512.5682 


-4.2696 


0.0125 


0.0250 


4512.6822 


-4.4391 


0.0107 


0.0530 


4512.7152 


-4.5081 


0.0103 


0.0580 


4515.6212 


-5.8950 


0.0099 


-0.0070 


4518.6525 


-4.5850 


0.0117 


0.0340 


4524.6750 


-5.9753 


0.0107 


0.0700 


4525.6755 


-4.0564 


0.0091 


0.0560 
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Figure 3. (a) The radial-velocity curve of WASP-14. A best- 
fitting MCMC model fit for the eccentric orbit is overplotted with 
the solid line, (b) Bisector spans as a function of RV. Residual 
about the MCMC fit as a function of transit phase (c), and as 
a function of time (d). The open triangle and filled circles are 
observations from the FIES and SOPHIE respectively. 



gorithm (see §3.3). The resulting normalized transit light 
curve of WASP-14 along with model fit is shown in Fig- 
ure [2jb) . From the model fit, we estimated a transit depth 
of about 12 mmag ((Rp/R*) 2 = 0.0102) and a rms residuals 
of 1.4 mmag around the best fit. 

2.3 Spectroscopic follow-up 

To derive the orbital parameters of the planetary system, 
we obtained spectroscopic observations with the Fibre-fed 
Echelle Spectrograph (FIES) mounted on the 2.5-m Nordic 
Optical telescope. A total of six high-resolution spectra of 
WASP-14 with an exposure time of 900 sec covering the 
wavelength region 4000-7350 A were obtained during 2007 
December 27-31 and 2008 January 25. A further 21 RV 
points were secured using the SOPHIE spectrograph on the 
1.93-m telescope at observatoire de Haute-Provence during 
2008 February 12-16, plus one point each on the nights 
of 2008 February 19, 22, 28 and 29. A detailed descrip- 
tion about the FIES a nd SOPHIE observations is given in 
l|Christian et al. 1 120081 ) . The journal of the FIES and SO- 
PHIE observations is given in Table[T] and includes barycen- 
tric Julian dates (BJD), RV measurements and associated 
error. A phase-folded radial- velocity curve overplotted with 
the maximum-likelihood orbit model is shown in the up- 
per panel of Figure [3ja) . This shows RV variations with a 



semi- amplitude K\ = 0.993 kms 1 and the period derived 
from the photometric observations. The velocity RMS to 
the orbital fit for FIES and SOPHIE are 5.8 and 11.0 ms" 1 
respectively. 

The variation of the line bisector with orbital phase 
is often used as a tool to dete ct blended stellar c ompan- 
ions or chromospheric activity l|Queloz et al.l [200 ll ). From 
the CCF we obtained the line bisector and calculated the 
bisector spans as noted in the last column of the Table [TJ 
and plotted, as a function of RV, in Figure |3Ib). We de- 
termined a signal-to- noise ratio SNR = 1.45 (see, Christian 
et al. 2008 for detail), indicating no significant correlation 
between the bisector span and RV variations. This supports 
the conclusion that the motion is indeed caused by a com- 
panion planet, which given the large RV semi-amplitude, 
suggests that WASP- 14b is a super-massive exoplanet. 

The RV residuals in Figure |3jc) show a small dispersion 
of o-(O-C) = 10.1 ms" 1 about the MCMC fit implying a 
X 2 = 23.62 for 27 degrees of freedom. This is consistent 
with the error bars on radial velocity measurements. We 
note this is also the order of magnitude of expected stellar 
'jitter' in the F5 dwarf host (see, e.g., Santos et al. 2000, 
Wright 2005). We have not noticed any systematic variation 
in the RV residuals during two months of our observations 
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Figure 4. The RV variations around the transit phase. Solid line 
shows best-fitting model fit, dashed line shows Keplerian model 
without Rositter-McLaughlin effect and dotted line shows the 
model fit of a transverse transit (A = 90°). 

(see, Figure |3jd)) hence do not suggest the presence of any 
third body in the system on the basis of present data. 



2.4 Rossiter-McLaughlin effect 

Several RV observations were secured with SOPHIE dur- 
ing a transit of WASP- 14b on the night of 2008 Febru- 
ary 13. The relatively high impact parameter of the system 
breaks the degeneracy between v sin i and the angle A be- 
tween the projected spin axes of rotation axis of the star 
and the orbit of the planet on the sky. We treated both 
these quantities as fitting parameters in the MCMC code, 
and verified that they were indeed uncorrelated. The result- 
ing vs'mi — 4.7 ± 1.5 km s _1 , determined from the appar- 
ent radial acceleration during the transit, is consistent with 
vsini = 4.9 ± 1.0 km s _1 derived from the spectroscopic 
analysis, and vsini = 3.0 ± 1.5 km s" 1 from the width of 
the CCF of the SO PHIE spectra which is calibrated follow- 
ing the method of ISantos et a l. (2002). We note here that 
FIES spectra and SOPHIE CCF provide a vsini averaged 
over all the stellar disk, whereas the R-M fit provides the 
v sini at the impact parameter of the transiting planet. The 
value of misalignment angle, A = —14^^ degrees, derived 
from the asymmetry of the R-M effect, is consistent with 
zero obliquity, but is of insufficient precision to rule out a 
substantial misalignment of the stellar spin and planetary 
orbital axes. 

An enlarged plot of RV variations around the transit 
phase is drawn in Figure U A continuous line is over-plotted 
for our best fit (wsini = 4.7 km s~ , A = —14°). We also 
overplot Keplerian fit by the dashed line determined from 
the RV measurements outside the transit in order not to 
be affected by the uncertainties of the R-M effect. One can 
see that the first four measurements performed in the first 
half of the transit are red-shifted by comparison to the Ke- 
plerian fit, whereas the last two measurements performed 
in the second half of the transit are blue-shifted by com- 
parison to the Keplerian fit. That ordinary feature for an 
aligned, prograde transit shows that the R-M anomaly is 
indeed significantly detected. A model of a transverse tran- 



sit drawn by dotted line in the same plot clearly exclude 
any variations due to A = 90°. A possible t ransverse tran- 
sit ha ve been detected in the case of XO-3b (|Hebrard et. al I 
l200ct). whereas HP 147 5 06b (HAT-P-2b) sho w aligned tran- 
sits (|Winn et al. 1120071 ; lLoeillet et al. 112008ft . This suggests 
that the history of WASP-14b is more similar to HAT-P-2b 
than the other massive planet XO-3b. Given the high eccen- 
tricity of the system, it would be interesting to make further, 
more densely-sampled observations of the R-M effect in this 
system. This will also allow to better constrain the system 
parameters, such as A. 

3 DETERMINATION OF SYSTEM 
PARAMETERS 

3.1 Tea from the IRFM 

We used archival TYCHO, DENIS and 2MASS magnitudes 
to estimate the host star's effe ctive temperature using th e 
Infrared Flux Method (IRFM) (|Blackwell fc Shallis ||1977m . 
IRFM yields a nearly model-independent determination of 
T e ff from the integrated stellar flux. We individually esti- 
mated T c ff for J, H and K bands and a weighted mean was 
determined for the final value of T e ff . In the case of WASP- 
14 we find T cff = 6480 ± 140 K. It should be noted that no 
extinction correction was applied in the IRFM fitting since 
the star is relatively nearby and out of the Galactic plane. 
The J — H — 0.22 and V — K = 1.13 colours also suggest 
similar estimat es for T c ff of ~ 6230 and 6440 K respectively 
l|Martin Il2006t ). This suggests the star is of spectral type of 
F5-7. We also determined an average angular diameter of 
9 = 0.075 ± 0.004 mas from the IRFM. 

3.2 Spectral analysis 

In order to perform a detailed spectroscopic analysis of the 
stellar atmospheric properties of WA SP- 14, we followed th e 
same methodology as described in IPollacco et al.l {2008). 
Our results from the spectral analysis are consistent with 
a F5 spectral type main-sequence star, in agreement with 
the F5-7 type determined from Infrared data, and the F5 
type listed in SIMBAD from the PPM catalog. In addi- 
tion we also analyzed a small region around the Li I 6708 
line to determine a lithium abundance that is defined as 
logiV(Li) = 12 + (Li/H) and estimated to be 2.84. The re- 
sulting parameters from our spectral analysis are given in 
Table [2] 

3.3 Model fit to determine planetary parameters 

Transit photometry combined with radial velocity measure- 
ments provides detailed information about the planetary or- 
bit and the stellar and planetary parameters. High precision 
photometric data was combined with the SOPHIE radial- 
velocity measurements in a simultaneous MCMC analysis to 
find the parameters of the WASP- 14 system and their covari- 
ance matrix. Our implementation o f MCMC for transiting 
exopla nets is presented in detail in ICollier Cameron et al.l 
1 20071). with an exten sion to fit eccentric orbits described in 
IPollacco et al.l l|2008h . Various parameters determined from 
model fit along with the corresponding 1-a error limits, are 
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Table 2. Stellar parameters for WASP-14. The last 6 parameters 
are derived from the spectral analysis of the FIES data. 



R.A. (2000) 


14 h 33 m 06.35 s 


DEC (2000) 


+21°53'40.9" 


V (mag) 


9.75 


B — V (mag) 


0.46±0.03 


Age (Gyr) 


~ 0.5-1.0 


Distance (pc) 


160±20 


Spectral Type 


F5V 


v sini (kms -1 ) 


4.9±1.0 


T«h (K) 


6475±100 


log g (cgs) 


4.07±0.2 


[M/H] (dex) 


0.0 ±0.2 


logiV(Li) 


2.84±0.05 
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Table 3. Planetary and stellar parameters for the WASP-14 sys- 
tem derived from the MCMC analysis along with lcr limits. 



Figure 5. Position of WASP-14 in R/M 1 / 3 - T eS plane. 
IGirardi et aTl (2000) solar metallicity stellar evolutionary tracks 
for different masses and isochrones of different ages are also 
shown. 



Parameter 



Symbol Values 



Transit epoch (HJD) 


To 


Orbital Period 


P 


Transit duration 


T 


Transit Depth 


(It) 2 


Impact parameter 


b 


Orbital Inclination 


i 


Orbital semi-major axis 


a 


Orbital eccentricity 


e 


Arg. periastron 


U! 


Stellar Reflex velocity 


K X 


Center-of-mass velocity 


7 


Stellar mass 


Af» 


Stellar radius 


R* 


Stellar density 


P* 


Stellar surface gravity 


logs* 


Planet mass 


M P 


Planet radius 


R P 


Planet density 


Pp 


Planetary surface gravity 


loggp 


Planet temp (A = 0) 


T P 


Photometric data points 


N df 


X 2 (photometric) 


~%phot 


Spectroscopic data points 


N df 


X 2 (spectroscopic) 


Xspec 



2454463.57583 ±0.00053 
2.243752 ±0.000010 days 
0.1275 

+0.6662 



days 



4 EVOLUTIONARY STATUS 

Stellar and planetary evolution 



T.l 



0.0102 
0.535 



+0.0028 , 

-0.0031 aa y s 

0.0003 ma S 
+0.031 r> 
-0.041 

+0.67 



84.32 deg 
0.036 ±0.001 AU 
0.091 ±0.003 
-106.629 ±o.678 rad 
0.993 ±0.003 km s" 1 
-4.990 ±0.002 km s" 1 



1 211 + ' 127 
-0.122 

1 306 + ' 066 
±.ouu _o.073 

542 + 079 

4 287 + ' 043 
-0.038 



M 
l: • 

P© 
[cgs] 



7 o 41 +0.508 
I -oil -0.496 
, 9S1 +0.075 
l.zoi _ .082 
o cm +0.636 
O.OU1 _ . 4 g 5 



Mj 

K.T 

, ,„ 

4.010 ™ [cgs] 
1866.12 tlllt K 

6482 

6452.78 

27 

23.62 



listed in Table E] Since a zero eccentricity orbit is expected 
for such a short-period planet system due to relatively small 
tidal circularisation timescale (see §4.2), we have also anal- 
ysed the model assuming a zero-eccentricity orbit which has 
increased x 2 by about 370. Using the stellar radius given in 
the table and angular diameter of the star determined from 
the IRFM, we estimated a distance of about 160 ± 20 pc for 
the host star. 



We estimated the age of WASP-14 using a maximum- 
likelihood fit to the stellar evoluti onary tracks for low- and 
intermediate-mass stars given by IGirardi et al~l l|2000l ). In 
Figured] we have shown position of WASP-14 in the R/M 1 ^ 3 
versus T e g plane. The evolutionary tracks of different stellar 
mass along with the isochrones of different ages taken from 
IGirardi et al~l (|2000h models are also drawn. This shows the 
star with T eff = 6475±100 K and a radius of 1.30±0.07 R , 
is consistent with being a main-sequence object and has an 
interpolated mass of 1.21 M Q from the models. The MCMC 
fit to the photometry and radial velocities (Table [3)l was 
computed using this value for the stellar mass, ensuring con- 
sistency between the evolutionary status and the system pa- 
rameters obtained from the MCMC fit. The resulting value 
of log g — 4.33±0.06 is higher but more precisely determined 
than that obtained from the spectral fit, and is consistent 
within the uncertainty of the latter measurement. Although 
there is significant uncertainty in the effective temperature 
of the star, WASP-14 seems to be consistent with being an 
approximate age of 0.5-1.0 Gyr. 

In F-type stars the convection zone does not penetrate 
deep enough to transport lithium to regions with tempera- 
tures above 6000 K, reducing the efficiency of lithium de- 
pletio n, with the exception of the Li gap or 'Boesgaar d 
—gap' (|Boesgaard k, Tripicco I Il986l ; iBalachandran I 119951 ). 
The temperature of WASP-14, T eff = 6475±100 K, places 
it close to the Li gap, where lithi um depletion is enha nced 
by an as-yet unknown mechanism ijBohm-Vitense 120041 ) . Al- 
though the rate of lithium depletion in the Li gap is not very 
well calibrated, the relatively high lithium content log AT (Li) 
= 2.84±0.05 and rotation speed vsxsxi — 4.9 ± 1.0 km s -1 
possibly indicate that WASP-14 is a young star. Similar 
high levels of lithium is also found in the Hyades cluster 
(age ~ 600 Myr) wit h log iV(Li) = 2.77±0.21 for stars with 
T cft = 6200 ± 150 K jSestito fc Randich 1120051 ). 

We further com pare our results with the models of 
iFortnev et al~l (|2007l ) where masses and radii are given for 
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a range of planetary masses, orbital separations and stellar 
ages. For WASP- 14b, these models predict a planetary ra- 
dius < 1.20 Rj for a host star of age 1 Gyr, while for an age 
of 300 Myr, the radius is predicted to be about 1.20-1.26 
Rj. A radius of 1.28 ± 0.08 Rj for the planet suggests that 
WASP- 14 possibly lies somewhere in the middle of these 
two age limits - in agreement with our earlier finding. For 
massive planets (> 5 Mj), the model shows <3% change 
in planet radius with cores of m ass from 0-100 Mq. In the 
context of iFortnev et al.l l|2007l ) model our observations do 
not place strong constraints on the core mass of WASP- 14b. 

4.2 Tidal evolution 

One of the interesting features of WASP-14b is its high or- 
bital eccentricity, e = 0.091, for its small orbital distance of 
0.036 AU. At this distance, any orbital eccentricity should 
induce strong tidal energy dissipati on within the planet lead - 
ing to circularisation of the orbit. Ma tsumura et all (|2008l ) 
determined a median eccentricity of 0.013 for close-in plan- 
ets with semi-major axis a < 0.1 AU. The nearly-circular 
orbits of most close-orbiting exoplanets indicate tidal cir- 
cularisation timescales significantly shorter than the system 
age. The high eccentricity of WASP-14b may thus indicate 
either a system age comparable to the tidal circularisation 
timescale, or the presence of an additional perturbing body 
in the system. 

In the tidal ev o lution scenario originally developed by 
Goldreich fc Soter I <l 19661) and more recen t ly dis cussed by 
Dobbs-Dixon et al. I (20041 '); IJackson et~al~l (120081 ) and oth- 
ers, the circularisation timescale (r e ) and spiral-in timescale 
(r a ) for a single planet are given by the logarithmic time 
derivatives of the orbital eccentricity and semi-major axis: 

T e = e/e; r a = a/a 

Using the equations given in IJackson et al. I (|2008l ) for 
close-in extrasolar planets, and the planetary and stel- 
lar parameters for the WASP-14 system, we obtain 
r e ~ l/[28.5/(Q p /10 5 ) + 6.6/(Q s /10 6 )] Gyr and r a ~ 
l/[0.5165(e/.089) 2 /(<2 P /i0 5 ) + 2.359/(Q s /10 6 )] Gyr. Here 
the tidal dissipation parameters Q v for the planet and Q 3 
for the star are assigned fiducial values of 10 5 and 10 6 re- 
spectively in acco rdance with the best-fi tting values derived 
from the study of IJackson et al.l (|2008l ). It should be men- 
tioned here that there could be considerable uncertainties 
in both these parameters; moreover, the timescales are in- 
versely proportional to the fifth powers of the planetary and 
stellar radii, both of which contain significant uncertainties 
in their determinations. 

From the expression for r e above, the circularisation 
timescale of the planet's orbit will be less than 1 Gyr if 
either Q p < 3.7 x 10 6 or Q s < 6.6 x 10 6 . The timescale for 
decay of the orbital semi-major axis of the planet via torques 
arising from the tidal bulge raised on the star could also be 
less than 1 Gyr if Q s < 2.4 x 10 6 . These same torques might 
be expected to spin the star up on a timescale of the order 
of 1 Gyr, since the angular momentum loss in the wind of 
a F star like WASP-14 is quite weak. A dissipation factor 
Qs > 6.6 x 10 6 could increase the life expectancy of the 
planet to a value similar to the main-sequence lifetime of 
the F-type host. 

If Q p is low and the circularisation timescale is short, 



the most plausible alternative mechanism for maintaining 
the high eccentricity of WASP- 14b is secular interaction 
with an additional planet in the system (see, e.g., Adams 
& Laughlin 2006). From the arguments above, however, it 
is not clear whether an additional planet is needed to main- 
tain the observed eccentricity over the ~ 1 Gyr lifetime of 
the star. This is only needed if the rate at which tidal energy 
is dissipated in the planet is high, with Q p < 3 x 10 6 . The 
slow stellar spin suggests that magnetic braking is dominat- 
ing over torques arising from the tidal bulge raised on the 
star. This in turn suggests a low rate of dissipation in the 
star, with Q s > 10 8 or so. A further discussion of eccentric 
orbits in cl ose-in exoplanets and the ir constraints on Q can 
be found in iMatsumura et al.l l|2008l ). 

Two other recently announced planets WASP-lOb 
(e=0.059; Christian et al. 2008) and WASP-12b (e=0.064; 
Hebb et al. 2008 ) also have apparently s ignificant orbital ec- 
centricity. While IShen fc Turner I ||2008T ) noted that sparsely 
sampled radial velocity curves are often interpreted as hav- 
ing significant orbital eccentricity, in the case of the WASP 
planets the spectroscopic data is usually obtained over a 
small number of cycles and with good precision. Conse- 
quently, we expect the eccentricities derived for these sys- 
tems to be free from bias or selection effects. 



5 DISCUSSION 

We have discovered a new massive exoplanet WASP-14b or- 
biting with a period of about 2.244 days around its host 
main-sequence star GSC 01482-00882. Spectral analysis of 
the star implies a spectral type of F5V with solar metallicity. 
High precision photometric and spectroscopic follow-up ob- 
servations reveal that the planet has a mass of 7.3±0.5 Mj, 
radius of 1.28±0.08 Rj suggesting a mean density of 3.5±0.6 
Pi (~ 4.6 gem -3 ), and an eccentricity of 0.091±0.003. The 
mean density of WASP-14b is high in compari son with a typ- 
ical H ot Jupiter density of 0.34-1.34 gem -3 IjLoeillet et al. I 
2008), and similar to that of rocky planets and makes it 
densest transiting exoplanets so far discovered with < 3 d 
orbital period. The plane t is too mass i ve to fit the mass- 
period relation given by iTorres et al. (2008), but its ra- 
dius iis_£OirsjsteirtwiHi_the theoretical radius expected from 
the IFortnev et al. ] (|2007l ) model. Spectral analysis reveals 
that WASP-14 has a high lithium abundance, logiV(Li) = 
2.84. This is consistent with its stellar temperature (T e a = 
6475±100 K) and age of about 0.5-1.0 Gyr which is further 
supported by the fitt ing of stellar evolutionary models of the 
iGirardi et aL~l(|20()Ch . 

High orbital eccentricities in close-in planets are often 
explained by perturbations from a third body (Jackson et al. 
2008 and references therein) or internal structures that re- 
sult in their tidal dissipation factor being significantly larger 
that that commonly assumed Q p ~ 10 6 (Matsumura et al. 
2008, Hebb et al. 2008). Although we do not see any sys- 
tematic variation in RV residuals during our two months 
of spectroscopic observations of WASP-14, we expect that 
long term radial velocity monitoring will help constrain the 
nature of the third body in this systems hence make it an im- 
portant target both for future transit-timing variation stud- 
ies and for longer-term RV monitoring to establish the mass 
and period of the putative outer planet. 
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WASP- 14b is one of the most mass ive transit- 
ing p lanets known alo n g with HAT-P-2b (|Bakos et al.l 
120071: IWinn et al. I [200^ ; lLoeillet et al. 1 l2008h and XO-3b 
Jjohns-Krull et al.ll20o3 IWinn et al. II2008I ; iHebrard et. al I 
2008) and its physical characteristics closely resemble with 
those of HAT-P-2 except the latter has a much longer or- 
bital period and smaller radius. However, there is still no 
firm explanation about the formation of such highly mas- 
sive planets (e.g., Hebrard et. al 2008). A quite interesting 
feature is that all these three massive exoplanets has un- 
usually large eccentric orbit for their short orbital period. 
Matsumura et al. (2008) argues that the new class of eccen- 
tric, short period transiting planets are still in the process 
of circularization and speculates that Q p could be as large 
as 10 9 for these planets. 

The success of theoretical models of planetary structure 
depends heavily on our precise knowledge of the basic phys- 
ical parameters of the planetary systems. Transit surveys in 
recent times have produced a large sample of transiting plan- 
ets that show a remarkable diversity in their mass, radius 
and internal structure. While most of the Jupiter-mass tran- 
siting planets are well explained by existing models, planets 
which show excessive mass for their small radius or relatively 
low mass for their large radius are yet to be explained sat- 
isfactorily by any available model. WASP-14b is one of the 
few massive exoplanets, some of which are even in closer or- 
bits than most other hot Jupiters. It poses a great challenge 
for theoretical models to explain their internal structure, 
atmospheric dynamics and heat distribution. 
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